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Abstract. The classical theory of linear response applies to 
statistical mechanics close to equilibrium. Away from equi- 
librium, one may describe the microscopic time evolution by 
a general differentiable dynamical system, identify nonequilib- 
rium steady states (NESS), and study how these vary under 
perturbations of the dynamics. Remarkably, it turns out that 
for uniformly hyperbolic dynamical systems (those satisfying the 
"chaotic hypothesis"), the linear response away from equilibrium 
is very similar to the linear response close to equilibrium: the 
Kramers-Kronig dispersion relations hold, and the fluctuation- 
dispersion theorem survives in a modified form (which takes into 
account the oscillations around the " attractor" corresponding to 
the NESS). If the chaotic hypothesis does not hold, two new 
phenomena may arise. The first is a violation of linear response 
in the sense that the NESS does not depend differentiably on 
parameters (but this nondifferentiability may be hard to see ex- 
perimentally). The second phenomenon is a violation of the 
dispersion relations: the susceptibility has singularities in the 
upper half complex plane. These "acausal" singularities are ac- 
tually due to "energy nonconservation" : for a small periodic 
perturbation of the system, the amplitude of the linear response 
is arbitrarily large. This means that the NESS of the dynamical 
system under study is not "inert" but can give energy to the 
outside world. An "active" NESS of this sort is very different 
from an equilibrium state, and it would be interesting to see 
what happens for active states to the Gallavotti-Cohen fluctua- 
tion theorem. 
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0. Introduction. 



The purpose of the present paper is to review the mathematics of hnear response in 
the framework of the theory of differentiable dynamical systems. 

Linear response theory deals with the way a physical system reacts to a small change 
in the applied forces or the control parameters. The system starts in an equilibrium or a 
steady state p, and is subjected to a small perturbation X, which may depend on time. In 
first approximation, the change 5p of p is assumed to be linear in the perturbation X. 

Apart from the linearity of the response X ^ 5p, one can make various physical 
assumptions: time translation invariance, time reversibility, causality (the cause precedes 
the effect), energy conservation, closeness to equilibrium. One can also find relations 
between the response to external perturbations and the spontaneous fluctuations of the 
system. Studying the consequences of the above assumptions for bulk matter has yielded 
the Onsager reciprocity relations, the Kramers-Kronig dispersion relations, the Green- 
Kubo formula and the fluctuation-dissipation theorem. Note that similar ideas have been 
used in the study of electrical circuits, in optics, and in particle scattering theory. It is 
also possible to discuss higher order, i.e., nonlinear response. 

If the physical system in which we are interested is described by classical mechanics 
(with external forces, and a deterministic thermostat [17], [20]), its (microscopic) time 
evolution is given by an equation 

| = ^W (0.1) 

in phase space M. We want to discuss the corresponding mathematical situation of a 
smooth (=differentiable) dynamical system (/*) on a compact manifold M. In the case of 
continuous time, (/*) is called a flow and is determined by (0.1) and x{t) = /*x(0), but we 
shall also consider the case of discrete time, where /"^ is the n-th iterate of a differentiable 
map f : M M. 

The study of linear response for general smooth dynamical systems (/*) encounters a 
number of difficulties, and we shall obtain both positive and negative results. In Section 1 
we discuss how physical notions (like equilibrium, entropy production, etc.) and principles 
(like causality, energy conservation, etc.) can be related to mathematical concepts per- 
taining to smooth dynamics. Then we shall analyze linear response for smooth dynamical 
systems in a number of different situations, both informally (Section 2), and rigorously 
(Sections 3, 4). 

If we have quantum systems instead of classical systems as considered here, the theory 
of nonequilibrium is in part similar and in part very different (one cannot use finite systems, 
and one looses the smooth dynamics on a compact manifold). We refer to V. Jaksic and 
coworkers (work in progress) for a comparison of classical and quantum nonequilibrium. 
Note that there is a vast literature on linear response that we have not quoted. Relevant 
to the approach discussed here is work by J. -P. Eckmann, C.-A. Pillet, G. Gallavotti, J.L. 
Lebowitz, H. Spohn, D.J. Evans, G.P. Morriss, W.G. Hoover, among others. 

Acknowledgments. 
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The work reported here was done in part at the ESI in Vienna, and in part at the 
Math Dept of McGill, Montreal. 

1. Smooth dynamics and physical interpretation. 

Let (/*) be a smooth dynamical system on the compact manifold M. Interpreting (/*) 
as time evolution, we describe a physical state by an (/*) -invariant probability measure 
p, and we assume that p is ergodic*. The invariance of p expresses the time translation 
invariance of our physical system. Certain simple time evolutions (completely integrable, 
for instance) turn out to be pathological from the physical point of view that interests us 
[48] , and it is necessary to make some chaoticity assumption saying that the time evolution 
is sufficiently complicated to avoid the pathologies. Let us be more specific. One can show 
that an infinitesimal change dx{0) in initial condition gives a later change dx{t) ~ exp Xt 
where A is called a Lyapunov exponent (if dimM = d, there are d Lyapunov exponents 
associated with an ergodic measure p). A weak chaoticity assumption is that there is 
at least one Lyapunov exponent A > 0. A strong assumption of this sort is the chaotic 
hypothesis of Gallavoti-Cohen [19], which says that (/*) is uniformly hyperbolic. This 
will be explained in Section 3 (but we shall also consider systems that are not uniformly 
hyperbolic). 

For general (/*) there is usually no invariant measure p absolutely continuous with 
respect to Lebesgue on M {i.e., such that p has density p{x) with respect to the Lebesgue 
measure dx in local charts of M). If there is an invariant measure p smoothly equivalent** 
to Lebesgue in local charts, and if it is ergodic, we say that p is an equilibrium state 
(this generalizes the situation where M is an "energy shell" H{p,q) — constant for some 
Hamiltonian H, and p is the corresponding normalized Liouville measure on M, assumed 
to be ergodic; p is also known as the microcanonical ensemble). A chaotic dynamical 
system (/*) typically has uncountably many ergodic measures. Which one should one 
choose to describe a physical system? A physically observed invariant state is known 
as a nonequilibrium steady state (NESS), and one can argue that it can be identified 
mathematically as an SRB probability measure, or SRB state. The SRB states have been 
defined first in the uniformly hyperbolic case [46], [35], [10], and then in general [25], [26]. 
At this point we do not give a formal definition but state a consequence (which holds under 
some extra condition***): if the ergodic measure p is SRB, there is a probability measure 

* i.e., there is no nontrivial invariant decomposition p = api + (1 — a)p2. Suppos- 
ing that the p-integral of continuous functions A is given by time averages: p{A) = 
limy^oo ^ Jq dt A{f^xo), then p is almost certainly ergodic by Bogolyubov-Krylov the- 
ory (see for instance Jacobs [21] Section 11.3). 

** This means that the density p{x), and l/p{x) are differentiable functions of x in local 
charts. 

*** The extra condition is, in the discrete time case of a system generated by a diffeomor- 
phism /, that p has no vanishing Lyapunov exponent, in the continuous time case (fiow) 
that there is only one vanishing Lyapunov exponent, corresponding to the direction of the 
fiow. See L.-S. Young [51] for further discussion. 
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i absolutely continuous with respect to Lebesgue (but in general not invariant) such that 



p{A) = / p{dx) A{x) = lim 




oo T 



1 



[ dt [ i{d x) A{f x) 



(1.1) 



when ^ is a complex continuous function on M (A is physically interpreted as an observ- 
able). This says that an SRB measure p is obtained from Lebesgue measure by a time 
average when the time tends to +00. The choice of +cxo (not —00) introduces a time 
asymmetry which will turn out to play in the mathematical theory the role played by 
causality as a physical principle. Note that equilibrium states are SRB states, and general 
SRB states come as close to equilibrium states as is possible when there are no absolutely 
continuous invariant probability measures. 

If dx denotes the volume element for some Riemann metric on M (or if dx is smoothly 
equivalent to Lebesgue in local charts), we define the entropy of an absolutely continuous 
probability measure £{dx) = £{x)dx by 



[In the formalism of equilibrium statistical mechanics, where dx is the Liouville volume 
element, S{£) is the Gibbs entropy associated with the density £{■)]■ Define such that 
{f*i){A) = i{A o /*) and write { f*i){dx) = it{x)dx, then S{£t) depends on t, and in the 
case of the time evolution (0.1) one finds 



where the divergence is taken with respect to the volume element dx. One can argue that 
minus the above quantity is the rate at which our system gives entropy to the rest of the 
universe. In other words, the entropy production e{£t) by our system is the expectation 
value of — divA", i.e., the rate of volume contraction in M. This definition extends to 
probability measures that are not absolutely continuous. In particular, (1.1) shows that 
e(p) = p(— divA:") is the appropriate definition of the entropy production in the SRB state 
p. Note that, since p is invariant, e{p) does not depend on the choice of the volume 
element dx (contrary to S{£)). The identification of the entropy production rate e(p) (a 
physical quantity) with the expectation value of the phase space volume contraction rate (a 
mathematical quantity) is important*. It permits in particular a study of the fluctuations 
of the entropy production, leading to the fluctuation theorem of Gallavotti and Cohen [19]. 

We perturb the time evolution (/*) in the continuous time case by writing 





d 
It 




dx 



X{x)^Xt{x) 



* The first reference we have for this identification is Audrey [1] 
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instead of (0.1). In the discrete time case, the smooth map / is replaced by / + o / 
(i.e., fx is replaced by fx + Xt{fx), t integer). The perturbation Xi is assumed to be 
infinitesimal, and may depend on the time t. 

Wc discuss first the continuoiis time case, assuming that a linear response 5tp is 
defined, and proceeding with the usual physical arguments (see Toll [49]). We take the 
expectation value 5tp{A) of an observable A and assume for simplicity that Xt{x) = 
X{x)(f){t). Linearity and time translation invariance then imply the existence of a response 
function k such that 

Stp{A) = J dt' K{t - t')^{t') 

The Fourier transform of 5tp{A) is then 

j dte''^*5tp{A) = k{uj)^{ijo) 

where the Fourier transform k of the response function is called the susceptibility. Note 
that, since the right-hand side is a product, there are no frequencies in the linear response 
that are not present in the signal 0. (Nonlinear response, by contrast, introduces harmonics 
and other linear combinations of the frequencies present in the signal). If (/> is square- 
integrable, J duj \(f){uj)\'^ may, in many physical situations, be interpreted as the energy 
contained in the ingoing signal. If we assume that our system does not increase the energy 
in the signal {conservation of energy) we see that the susceptibility must be bounded: 
\k{-)\ < constant. Note that our physical assumption of "energy conservation" need not 
apply to a general dynamical system. In the discrete time case, the situation is analogous 
to that just described, Fourier transforms are replaced by Fourier series and it is convenient 
to introduce the variable A = e*^ so that the susceptibility is replaced by a function ^'(A). 

In our physical discussion, causality is expressed by the fact that K,{t — t') vanishes 
when t < t' . This, together with the boundedness of the susceptibility, implies that k{-) 
extends to a bounded analytic function in the upper half complex plane, and that the real 
and imaginary parts of k{uj) (for uj £ R) satisfy integral relations known as the Kramers- 
Kronig dispersion relations (see the discussion in [49]). When dealing with a general 
dynamical system, causality is replaced by the assumption that the state p of our system 
is an SRB state. To clarify this point we shall in the next section make a nonrigorous 
calculation of 5tp, and see how causality appears. We shall also discuss the special case 
when p is an equilibrium state (perturbation close to equilibrium) and understand how the 
fluctuation- dissipation relation arises. Away from equilibrium only part of the fluctuation- 
dissipation relation will survive. 

2. Linear response: an informal discussion. 

We shall now evaluate the linear response 5p by a nonrigorous calculation (using what 
is called first order perturbation theory in physics). Le us consider a discrete dynamical 
system (f^) where / : M ^ M is a smooth map, and the formula (1.1) is replaced by 

p{A) = lim - V / e(dx)A{f''x) = lim - V / (/'=*£) (da;) A(x) (2.1) 

n—^oo n / n— >oo n I 
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In (2.1), we assume that £{dx) = i{x)dx is an absolutely continuous probability measure 
on M. The formula (2.1) holds if p is an SRB measure for a diffeomorphism / of M, and 
also if p is an absolutely continuous invariant measure (a.c.i.m.) for a map / of an interval 
[a, b] C R. Replacing fx by fx — fx + X{fx) we have to first order in X 

k 

~fx = fx + Y,{Tpj'-^)X{rx) 

j=l 

where T^f denotes the tangent map to / at the point x. Thus 

k 

ACf^x) = Aifx) + A'ifx) J2iTpxf-')Xifx) 

k 

= A{fx) + X{fx) ■ Vp,{A o f-^) 



hence 



n k 



1 " /• 

6p{A)^ lim h{dx)X{fx)-Vp,{Aof-^) 



k=i j=i 

n k 



1 /■ 

= hm UfH{dx))X{x)-V,{Aof^-^) 



k=l j=l 



= [{ri{dx))X{x)-W,{Aof) 

n—^oo n I 
i>0 j=l 

If we interchange in the right-hand side lim^^oo ^ind X]i>o (without a good mathematical 
justification!), and use lim„_,oo ^ Y^^Zl f^*^ = p we obtain formally 

oo „ oo „ 

5p{A)^Y. / p{dx)x{x)-v,{Aor) = Y, / p{dx)x{r''x)-Vf-r.Mon (2.2) 

n=0 n=0 

The physical meaning of this formula is that the change of p{A) due to the perturbation 
X is a sum over n of terms corresponding to the perturbation acting at time —n. The 
fact that the sum extends over n > may be interpreted as causality, and results from the 
asymmetry in time of the formula (2.1) defining the SRB measure p. Replacing X{f~'^x) 
by e^"'^X{f~"'x) in the right-hand side of (2.2) we obtain the susceptibility which, as a 
function of A = e*"''*', is 

oo „ 

^(A) = ^ A" / p{dx) X{x) ■ V,(A o D (2.3) 
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Since / has bounded derivatives on M, the power series in A defined by the right-hand 
side of (2.3) has nonzero radius of convergence. FormaUy, Sp{A) = ^'(l), but we are not 
assured that \E'(1) makes sense. 

In the continuous time case (2.3) is replaced by 

k{iu) ^ dt e"^* j p{dx) X{x) • V^(A o f) (2.4) 

and formaUy 5p{A) — k(0) (this corresponds to taking Xt{x) = X{x) or (f)(t) = 1 in 
Section 1) but we are not assured that k(0) makes sense. Comparing with the physical 
discussion of Section 1, we note that K,{t) = for t < 0, i.e., causality is satisfied. But 
it may happen (sec Section 4) that K{t) = f p{dx) X (x) ■ VxiA o /*) grows exponentially 
with t, so that k{uj) does not extend analytically to the upper half plane. This apparent 
"violation of causality" in fact means that "conservation of energy" is violated: when hit 
by the periodic perturbation e'^^^X, the system may give out (much) more energy than it 
receives. 

We consider now the situation where p is an equilibrium state. Thus p{dx) — p{x)dx 
in a local chart and we may define diVpX by 



It is convenient to write simply dx for the volume element p{x)dx and div^jX for diVpX^x). 
We obtain then 

J p{dx) X{x) .Vx{Aor) = - j dx (div,X)yl(/"a;) 

Notice that the right-hand side is a correlation function in the time variable n. If we 
assume that this correlation function tends exponentially* to when n — > oo it follows 
that the radius of convergence of ^'(A) is > 1, and (2.3) becomes 

oo » 

*(^) = - 5^ AM da; (div,X)^(ra;) 

n=0 

which makes sense for A = 1. Similarly we obtain from (2.4) in the continuous time case: 

k{uj) = - dt e*'^* J dx (div^X)yl(/*a;) (2.5) 

The susceptibility k{Lv) appearing in the left-hand side of (2.5) gives the linear response of 
our dynamical system to a periodic signal which puts the system outside of equilibrium, 



* Note that J p{dx) diVy^X (x) = so that the correlation function tends to when 
n — > oo if the time evolution is mixing. 
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i.e., in a so-called dissipative regime. The right-hand side is constructed from a time 
correlation function which describes the fluctuations of our system in the equilibrium state. 
The relation between dissipation and fluctuations expressed by (2.5) is a form of the so- 
called fluctuation- dissipation theorem ((2.5) is also related to the Green-Kubo formula). 
A physical interpretation of the fluctuation-dissipation theorem is that kicking the system 
outside of equilibrium by the perturbation X is equivalent to waiting for a spontaneous 
fluctuation that has the same effect as the kick. The reason that this is possible is the 
absolute continuity of p. 

If we assume that the correlation function in the right-hand side tends exponentially* 
to when t ^ oo, we see that k{-) extends to an analytic function in {cu G C : Imu > — e} 
for some e > and that the real and imaginary parts of k{-) on R are related by Hilbert 
transforms: these relations are the Kramers-Kronig dispersion relations (see Toll [49]). 

We return now to the study of the susceptibility (2.3) or (2.4) when p is an SRB but 
not necessarily an equilibrium state. At this point we need a brief description of the ergodic 
theory of smooth dynamical systems following the ideas of Oseledec (see [31], [36]), Pesin 
[33], [34], Ledrappier, Strelcyn, and Young (see [25], [26]). For deflniteness we discuss the 
discrete time case of a diffeomorphism / : M — > M where M has dimension d. Given an 
ergodic measure p for /, there are d Lyapunov exponents Ai < . . . < which give the 
possible rates of exponential separation of nearby orbits (almost everywhere with respect 
to p). For p-almost every x, there are a stable (or contracting) smooth manifold V| and 
an unstable (or expanding) smooth manifold V" through x. The dimension of is the 
number of Lyapunov exponents < 0, and the dimension of is the number of Lyapunov 
exponents > 0. The manifold V| is shrunk exponentially under iterates of / while is 
shrunk exponentially under iterates of f~^. The ergodic measure p is SRB if and only 
if it is smooth along unstable directions. This may be taken as a general definition of an 
SRB measure, and means that there is a set S with p{S) = 1, and a partition of S into 
pieces Eq, C of unstable manifolds such that the conditional measure cTq, of p on 
is absolutely continuous with respect to Lebesgue on V^. The manifolds V|,V^ do not 
depend continuously on x (only measurably). If and Vy are unstable manifolds with y 
close to X, one can define a holonomy map tt of part of Vy to part of along the stable 
manifolds (a stable manifolds through a point of Vy may "curve back" before it hits 
but, in terms of Lebesgue measure, most of Vy is mapped to most of if y is sufficiently 
close to x), and the map tt is absolutely continuous**: sets of Lebesgue measure are sent 
to sets of Lebesgue measure 0. 

* Exponential decay of correlations does not always hold for smooth dynamical systems, 
even if they are uniformly hyperbolic, but is still a natural assumption [15]. In the statisti- 
cal mechanics of bulk matter, time correlation functions decay more slowly, say like 
which holds for diffusion in u dimensions, and it is natural to assume absolute integrability 
in time, so that the susceptibility is analytic only in the upper half plane. Incidentally, 
this means that bulk matter with dimension < 3 is expected to behave pathologically 
with respect to linear response. 
** This is a form of Pesin's theorem of absolute continuity of foliations. See for instance 
[6] p.302. 
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Suppose now that the SRB measure p for / has no vanishing Lyapunov exponent. 
We can then write (for p-almost all x) X{x) = X^{x) + X^{x) where X^[x) is in the 
stable direction (tangent to V|) while X'^[x) is in the unstable direction (tangent to V^). 
Inserting this in (2.3) we find 

oo „ oo „ 

*(A) = 5^ A- / p{dx)X\x) • V,(A o r) + 5^ A- / p{dx)X^{x) • V,(A o /-) 

n=0 n=Q 

The X*-integral may be rewritten as 



j p[dx)[(Tf-r.J^)X\r^x))-V,A 



where {Tf^)X^ decreases exponentially with n. Since p is an SRB measure, the X^- 
integral may be rewritten in terms of integrals with respect to aa{dx) on pieces Eq. of 
unstable manifolds V^, where aa{dx) is absolutely continuous with respect to Lebesgue on 
V^. Introducing a divergence div" in the unstable direction one may hope to rewrite the 
X"-integral as 

' p{dx){dw:^x-{x))A{rx) 



This is a correlation function with respect to the time variable n, and one may hope that it 
tends to when n ^ oo. In conclusion one may hope that ^'(A) has radius of convergence 
> 1 and that 6p{A) = *(1). 

In the continuous time case, the hope is that we may rewrite (2.4) as 

k%co) = dt e-* j p{dx){{Tf-,j')X%f-'x)) ■ V^A 

k'^iu) = -J^^dte^^'J p{dx){diyZX^''{^))A{fx) 

where X^'^ is the component of X in the center-unstable direction corresponding to the 
Lyapunov exponents > (one zero exponent for the "flow direction", i.e., the direction 
of X). An optimistic guess would be that k{iv) extends to a holomorphic function for 
Imu; > — e with e > 0, and 6p{A) = k{0). Note that k'^^ is formally the Fourier transform 
of a time correlation function (cut to t > 0), i.e., k^"^ formally conforms to the fluctuation- 
dissipation theorem (as in the case where p is an equilibrium state). In order to interpret 
remember that, if p is not an equilibrium state, p is singular with respect to Lebesgue. 
One might say that p is concentrated on an attractor* ^ M. A perturbation that kicks the 



* The idea that p is concentrated on an attractor C M is geometrically appealing, but 
in fact the support of p may be the whole of M. The geometric notion of an attractor in 
M should thus be replaced by the idea of p as a measure-theoretic attracting point. 
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system in the stable direction away from the attractor is not equivalent to a spontaneous 
fluctuation. The effect of such a perturbation (the system oscillates and tends to the 
attractor) is described by [39], [28]. 

In Section 3 we shall see that in the uniformly hyperbolic case things work out basically 
as indicated above. In general however, p{A) is not expected to depend differentiably 
on the dynamical system. In order to understand what happens to linear response in 
nonhyperbolic situations we shall discuss in Section 4 the case of unimodal maps / of the 
interval and see how p{A) may depend nondifferentiably on /, yet look differentiable for 
the purposes of physics. 

3. Linear response: the uniformly hyperbolic case. 

In this Section we shall discuss the case of a dynamical system (/*) restricted to a 
neighborhood U of a, hyperbolic attractor* K C M. (This includes the situation that 
/ is an Anosov diffeomorphism of M or (/*) an Anosov flow on M, in those situations 
K = U = M). For completeness we now give a certain number of deflnitions. These 
definitions make use of a Riemann metric on M, but do not depend on the choice of the 
metric. 

If / is a diffeomorphism, we say that i^T is a hyperbolic set for f if TkM (the tangent 
bundle restricted to K) has a continuous r/-invariant splitting TkM = E'^ (B and 
there are constants c, A > such that 

||T/"v|| < ce-"^||v|| if veE',n>0 

< ce-'^^||^;|| if veE'',n>0 
{E^ and E'^ are called the stable and unstable subbundles of TkM). 

If (/*) is the flow associated with a vector fleld X (see (0.1)), we assume (for simplicity) 

that X does not vanish on K and we say that i^T is a hyperbolic set for (/*) if TkM has 
a continuous invariant splitting TkM = E"^ (B E^ (B E^, and there are constants c, A > 
such that 

E'^ is one-dimensional and E^ = 'RX{x) 

||T/*z;|| < ce-^*||^;|| if veE',t>0 
llTZ-^-yll < ce-^^ll-ull if veE'',t>0 

We say that K is a basic hyperbolic set (for a diffeomorphism or flow) if 

(a) K is a compact invariant set as above, 

(b) the periodic orbits of f*\K are dense in K, 

(c) is topologically transitive, i.e., the orbit {f*x) of some x E K is dense in K, 

* Hyperbolic attractors are also called Axiom A attractors [47] . Hyperbolicity as deflned 
here is uniform hyperbolicity. Weaker forms of hyperbolicity are not considered in this 
Section. 
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(d) the open set U D K can be chosen such that Htf^U = K. 

In particular we say that is a hyperbolic attractor (or Axiom A attractor in the 
terminology of Smale [47]) if one can chose U such that f^UcU for all sufficiently large 
t. We have then K = r)t>ofU. 

Uniform hyperbolicity has been much studied in relation with structural stability, to 

a small perturbation of (/*) corresponds a small perturbation of the hyperbolic set K. In 
particular, if the diffeomorphism / has a hyperbolic attractor K, and / is close to / in a 
suitable topology, then / has a hyperbolic attractor K and there is a homeomorphism 
h : K ^ K close to the identity such that f\K = hofoh~^. In the case of a flow 
(/*) there is a similar result, but a reparametrization of f*\K is necessary. An important 
tool in the study of uniformly hyperbolic dynamical systems is constituted by Markov 
partitions (introduced by Sinai [44], [45], with important contributions by Bowen [7]. [8], 
[9]). Using Markov partitions, it is possible to replace problems about measures on a 
hyperbolic attractor by problems of equilibrium statistical mechanics of one-dimensional 
spin systems. The transition is via symbolic dynamics, and makes available some efficient 
tools like transfer operators (introduced by Ruelle, see in particular [2] which is a good 
source of references on earlier literature). The body of knowledge accumulated in this 
direction is known as thermodynamic formalism (see [37], [32]). 

Of interest to us here are the results concerning SRB measures. On a hyperbolic 
attractor K there is exactly one SRB measure p. Writing p{A) — J p{dx)A{x) we may 
characterize p as follows (with U D K as above). 

(i) For Lebesgue-almost every x E U, 



if A is a continuous function U ^ C. 

(ii) If £{x)dx is a probability measure with support in U and absolutely continuous 
with respect to Lebesgue, then 



in the diffeomorphism case, and similarly in the flow case. 

There is a rich geometric theory of hyperbolic sets (beginning with [47]) which we 
cannot properly describe here. Let us mention that for every point x of a hyperbolic set 
there are a stable and an unstable manifold V|, V^, exponentially contracted by /* or /~* 
when t — > oo. If K is a hyperbolic attractor, then V" C K. In the flow case, it is useful 
to consider the center- unstable manifold V^", union over t of the /*V^. The SRB measure 




N-l 



(diffeomorphism case) 




(flow case) 




7V-1 
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p is the only ergodic measure on K such that its conditional measures on unstable (or 
center-unstable) manifolds are absolutely continuous with respect to Lebesgue on those 
manifolds (see Section 2). 

There are other important characterizations of the SRB measure p that will, however, 
not be used here. One of them is that the invariant* fiKsip) (rate of creation of information 
in the state p) is equal to the sum of the positive Lyapunov exponents of p (this sum is 
also equal to the expectation value in p of the logarithm of the unstable Jacobian = rate of 
growth of the unstable volume element). Another is that p is stable under small stochastic 
perturbations of the dynamics. 

We now have the concepts that allow us to discuss the dependence of the SRB measure 
p on the dynamical system (/*). We first discuss the discrete time case. 

Let Kq be a hyperbolic (i.e., Axiom A) attractor for the diSeomorphism fo, and 
suppose for simphcity that fo\Ko is mixing**. If f is allowed to vary in a small neighbor- 
hood of fo, there is a hyperbolic attractor K for f, depending continuously on f and a 
unique SRB measure p for f with support K. Furthermore 

(a) there is a neighborhood N of fo such that if A : M ^ H is C^, then f i— > p{A) 
is differentiable in M , 

(b) the first-order change Sp{A) when f is replaced by f -\- X o f is given by Sp{A) = 
^{l), where the power series 



has a radius of convergence > 1, 

(c) if X^jX"^ are the components of X in the stable and unstable directions, we may 
write \E'(A) = \E'''(A) + \E'"(A), where the power series 



both have radius of convergence > 1. 

The proof of (a) appears in [24], while (a),(b),(c) are proved in [38]. 

* See [5]. The Kolmogorov- Sinai invariant Hks is also known as entropy, but should 
not be confused with the Gibbs entropy discussed in Section 1. 

** This is equivalent to requiring that K is connected. In general K would consist of 
m connected components permuted by /, and one reduces to the mixing situation by 
considering restricted to one of the connected components. 




oo 



n=0 





oo 
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The divergence in the unstable direction div^X" that appears above can be shown to 
be a Holder continuous function. Therefore, the coefficients of the power series are the 
values of a correlation function n J p{x) A{f"'x) B (x) tending exponentially to 0. One 
can also show that \E'"(A) extends meromorphically to a circle |A| < -R with R > 1; its 
poles are Ruelle-PoUicott resonances (see [2]) corresponding to fluctuations of the system 
(if*), p) in accordance with the fluctuation-dissipation theorem. The poles of \E''^(A) would 
correspond to resonances in the "oscillations of the system around its attractor". Near 
equilibrium, i.e., when p is absolutely continuous and we may write p{dx) = dx, we find 
that the coefficients of are also the values of a correlation function, and the two kinds 
of poles become the same. It would be interesting to discuss examples where the "stable" 
resonances separate from the "unstable" resonances as one moves away from equilibrium. 

The following figure shows the singularities of ^'(A) in the complex A-plane: no pole 
for |A| < 1, some poles (crosses) for 1 < |A| < R, possible essential singularities for |A| > R. 




We discuss now the flow case (continuous time). 

Let the vector Geld X + aX on M deRne a Row (/*) with a hyperbohc (i.e., Axiom 
A) attractor Ka depending continuously on a & (— e,e). There is then a unique SRB 
measure pa for (/*) with support Ka- Furthermore 

(a) ifA-.M^R is C^, then a i-^ Pa{A) is on (-e, e), 

(h) the derivative dpa{A)/da is the limit when a; — > for u > of 

kaiu) = e'^'dt J Pa{dx)X{x) • V,(^ o fl) 

where ka{w) is holomorphic for Imco > 0, 

(c) the function u i— > ka{cij) extends meromorphically to {u : Ima; > —A}, and the 
extension has no pole at co = 0. 
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This is proved in [42], using the machinery of Markov partitions. A different proof 
has been given in [11] in the case of Anosov flows. 

Note that we did not require (/*) to be mixing, and ka{u!) may thus have poles on 

the real axis. If (/*) is mixing, ka{uj) might still have poles arbitrarily close to the real 
axis. To obtain a discussion analogous to that for diffeomorphisms, it is natural to define 
a center-unstable divergence C — div'^"(X'^ + X"), where we take a = for simplicity. If 
J dt\po{{A o /o)C)| < oo we have 

d f°° f 

^Pa(A)|a=0 = dtj po{dx)X{x)-V^{Aof^) 

If the correlation function t i— > po{{A o /o)C) tends to exponentially at oo, the poles 
of ko{u!) stay a finite distance away from the real axis. A number of results are known 
on the decay of correlation functions for hyperbolic flows (see in particular Chernov [13], 
Dolgopyat [14], [15], Liverani [27], Fields et al. [18]). 

The following figure shows the singularities of k{u}) in the complex cu-planc: no poles 
for Ima; > 0, some poles for > Ima; > —A, possible essential singularities for Ima; < —A 



i 

o 


^Imco 


X 

X 


> 

X Re 

X 


-A 





To compare with diffeomorphisms, remember that A = e*'^. 

Let us mention at this point a class of dynamical systems that are " almost" uniformly 
hyperbolic, namely the Lorenz system and related flows (see for instance Chapter 9 in [6]). 
Numerical studies [29] seem to indicate that the Lorenz system behaves like uniformly 
hyperbolic dynamical systems with respect to linear response. A mathematical study 
would here be very desirable. For partially hyperbolic systems see [16]. 



4. Linear response: the case of unimodal maps. 

Among dynamical systems that are not uniformly hyperbolic, the unimodal maps of 
the interval have been particularly well studied. These are smooth noninvertible maps 
f : I ^ I, where / = [a, b] is a compact interval in R. One assumes that a < c < b and 
that f'{x) > for a; e [a, c), f'{x) < for a; e (c, b]: this makes / one-humped (unimodal). 
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It is known since Jakobson [22] that many unimodal maps have an a.c.i.m. p{x)dx, i.e., an 
invariant measure absolutely continuous with respect to Lebesgue on [a, b]. (The concept of 
an SRB measure discussed earlier is replaced by that of an a.c.i.m. p{dx) = p{x)dx in the 
present 1-dimensional situation). Consider the specific example of maps fk ■ [0, 1] — > [0, 1] 
defined by fkX = kx{l — x). These have an a.c.i.m. pk for a set S of values of k with 
positive Lebesgue measure in (0, 4]. But one also knows that the complement of S is dense 
in (0,4]. How then could the function k ^ pk he differentiable? One idea [40] is to use 
differentiability in the sense of Whitney [50]: find a differentiable function (p such that 
(f){k) — Pk{A) when k belongs to some set S C 5 and define the derivative dpk{A)/dk to 
be (j)'(k) when k eT,. 

The above considerations suggest the following program: start with a unimodal map 
/ with an a.c.i.m. p, perturb f to f + X o f and define the corresponding derivative of 
p{A) to be 

6p{A) = 

where 

*(A) = J2^^ p{dx)X{x)-A{rx) 

n=0 

Then, show that 5p{A) is indeed a Whitney derivative in some sense. The conclusions of 
this program are as follows: 

(i) even with the idea of Whitney derivative, it appears that f ^ p is (mildly) non- 
difFerentiable, 

(a) the radius of convergence of \1/ is < 1, i.e., the susceptibility function oj i-^ ^'(e*'^) 
has singularities in the upper half-plane. 

We now give an idea of how such conclusions can be reached, referring to [41], [23], 
[43]* for the detailed assumptions and proofs. We note first that it is no loss of generality 
to suppose that the interval [a, b] and the critical point c satisfy fc — b, fb — a, and we 
assume a < fa < b. If the density p(-) of the a.c.i.m. is differentiable and nonzero at c, 
the invariance of p{x)dx under / implies that p{x) has a spike ~ (6 — near b. In 

fact, for each point f^b = f'^'^^c of the critical orbit (with n = 0, 1, . . .) there is a spike on 
one side of with singularity 

Cn\x-rb\-'/^ 

where 

^ n— 1 

Cn^p{c)\-nc)l[f'{fb)\-'/' 
k=0 



* Let us note that Baladi and Smania [3] have made a detailed study of linear response 
for piecewise expanding maps of the interval. The theory of these maps is in some respects 
very similar to that of smooth unimodal maps, in other respects quite different has no 
singularities inside the unit circle). 
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The easiest situation to analyze is when the critical point is preperiodic. Specifically 
we assume that f^c belongs to an unstable periodic orbit of period i for some /c, £. One can 
then prove that the radius of convergence of \I/(A) is < 1. In fact \E' has £ poles regularly 
spaced on a circle of radius < 1 (but no singularity at A = 1). The situation is thus 
analogous to that of a hyperbolic attractor, but with some poles inside the unit circle. 
This gives an easy example of our assertion (ii) above. 

We consider now a more general class of unimodal maps. First, we assume that / has 
a compact invariant hyperbolic set K. This is no big deal: we can assume that > 1 

for X away from some neighborhood of c, so that the points with orbits avoiding a suitable 
neighborhood of c automatically form a hyperbolic set K. Misiurewicz [30] has shown that 
if some point of the critical orbit is in K (say f^c G K), then / has an a.c.i.m. Here 
the singularities of inside the unit circle are expected to be worse then poles, probably 
forming a natural boundary. The interest of the Misiurewicz situation is that it can be 
perturbed: replacement of a Misiurewicz map /o (with hyperbolic set Kq, and critical point 
Co) by / replaces Kq by K, with a homeomorphism h : Kq — > K such that f oh — ho Jq on 
Kq. If f^c e K, then / is again a Misiurewicz map, with an a.c.i.m. p, and we can study the 
dependence of p on /. The evidence is that p does not depend differentiably on /. This is 
because the n-th spike has an amplitude decreasing exponentially like 1 11^=0 /'(Z'^^) 
but it moves around at a speed that can be estimated* to increase exponentially like 

Note that this evidence of nondifferentiability is not a proof! The nondifferentiability 
of / I— > p will probably not be very important in physical situations, because the spikes 
of high order (which lead to nondifferentiability) will be drowned in noise and therefore 
invisible. Using Misiurewicz maps we have thus argued [43] that, for such maps, / i— > p is 
mildly nondifferentiable, in agreement with our assertion (i) above. A study of the more 
general CoUet-Eckmann maps is under way [4]. 

To summarize, the study of unimodal maps reveals two new phenomena that should be 
present in more general dynamical systems: (i) nondifferentiability of / i— > p, and (ii) the 
apparently "acausal" singularities of the susceptibility in the upper half complex plane. 
Of these phenomena, (ii) may be most easy to observe (see for instance the numerical 
study by Cessac [12] of the Henon map). Both (i) and (ii) can occur only for dynamical 
systems that are not uniformly hyperbolic but, more specifically, they are related to "energy 
nonconservation" . This means that if the system is subjected to a periodic perturbation of 
small amplitude, the expectation value of observables may undergo a change of arbitrary 
large amplitude: the system is not passive or inert, it gives away energy to the outside 
world. To be specific, we can say that the dynamical system (/*) with SRB measure p is 
active if the corresponding susceptibility has singularities with Imu > 0. Such singularities 
are expected in systems where "folding" causes tangencies between stable and unstable 
manifolds, as happens for the Henon map. Indeed, the unimodal maps are a 1-dimensional 
model for folding in smooth dynamical systems. 

* This estimate fails if certain cancellations occur, and cancellations indeed occur when 
f^c = H/qCq. Those / such that f^c = H/qCq form the topological conjugacy class of /o, 
and the map f ^ p can be shown to be differentiable on a topological conjugacy class [43] . 
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5. Conclusions. 



To understand linear response for a NESS away from equilibrium, one is led to inves- 
tigating linear response for a general smooth dynamical system on a compact manifold. 
It turns out that for hyperbolic dynamical systems, the theory of linear response is sim- 
ilar to the classical theory close to equilibrium: The Kramers-Kronig relations hold, and 
the fluctuation-dissipation theorem is modifled by taking into account "oscillations of the 
system around its attractor". For nonhyperbolic systems, linear response may fail in the 
sense that the NESS does not depend differentiably on the parameters of the system. But 
this nondifferentiability may not be visible in physical situations. Perhaps more important 
is the failure of dispersion relations: there may be singularities of the susceptibility in the 
upper half complex plane. This happens when the system under consideration is active: 
physically this means that it can give away energy to the outside world. Clearly, the "close 
to equilibrium" paradigm has to be drastically revised in the case of active systems. It 
appears that both mathematical analysis and numerical simulations will be necessary to 
proceed to a new paradigm covering active systems, and to see in particular what happens 
to the Gallavotti-Cohen fluctuation theorem [19]. 
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